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A THEORETICAL ANALYSIS OF THE DYNAMIC LATERAL STABILITY 
AND CONTROL OF A PARAWING VEHICLE* 

By Joseph R. Chambers and Peter C. Boisseau 
Langley Research Center 

SUMMARY 

A theoretical analysis has been made to provide an understanding of some of the 
fundamentals of the dynamic lateral stability and control of parawing vehicles. The 
analysis is treated mainly in terms of the effects of vertical center-of-gravity position, 
since this factor is responsible for most of the unusual characteristics of such vehicles . 

Increasing the vertical distance between the payload and parawing led to decreased 
damping of the roll -subsidence mode and increased damping of the spiral mode. It also 
caused an increase in the damping and period of the lateral oscillation and a reduction in 
the ratio of roll to sideslip. The stability derivative C^p (yawing-moment coefficient 
due to rolling angular velocity) and the mass parameter Kxz (nondimensional product- 
of-inertia parameter) were found to be the principal factors influencing the damping of the 
lateral oscillation as the center-of-gravity location was lowered with respect to the 
parawing. 

There was little effect of vertical center-of-gravity location on the initial roll 
response to wing bank control but as the motion progressed, the configurations having the 
lower center-of-gravity locations displayed the greater roll response. The effects of 
vertical center-of-gravity location on the mass distribution of the system produced sig- 
nificant changes in the initial yawing motion following the application of wing bank control. 
It was found that the large nose -down inclination of the principal axis of least inertia 
caused by lowering the center-of-gravity position relative to the parawing tended to 
reduce the initial adverse yawing motion of the system. 

INTRODUCTION 

Recently, the National Aeronautics and Space Administration has been conducting 
general investigations to provide some fundamental information on the stability and con- 
trol characteristics of parawing configurations. (See refs. 1 to 4, for example.) The 
present investigation has been made to improve the basic understanding of the dynamic 
lateral stability and control of parawing configurations. The characteristics of such 

*The material presented herein includes information from a thesis entitled "An 
Investigation of the Dynamic Lateral Stability and Control of a Parawing Vehicle" by 
Joseph R. Chambers submitted to Virginia Polytechnic Institute, Blacksburg, Va., May 
1966 in partial fulfillment of the requirements for the degree of Master of Science. 
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configurations would be expected to differ from those of conventional aircraft in signifi- 
cant ways because certain physical characteristics of the configurations which are impor- 
tant to dynamic stability differ markedly from the characteristics of conventional air- 
craft. These physical characteristics are: (1) a center-of-gravity location far below the 
wing, (2) a mass distribution in which the mass is distributed mainly along the vertical 
axis, and (3) a wing having a shape which tends to give unusually large values of lateral 
force due to sideslip. The second factor is a consequence of the first, and the third fac- 
tor is important because of the first. Consequently, since the low center-of-gravity posi- 
tion was expected to be the principal cause of differences in stability of the parawing con- 
figuration from conventional aircraft, the investigation of stability is treated mainly in 
terms of the effect of vertical center-of-gravity location on lateral stability. The effects 
of individual stability derivatives were also investigated to determine those that were of 
predominant importance and to determine those that, in changing as the vertical location 
of the center of gravity was changed, were primarily responsible for the effects of 
center-of-gravity location on dynamic lateral stability. 

The investigation consisted of a theoretical determination of the dynamic lateral 
characteristics of an unpowered parawing-payload combination gliding at maximum lift- 
drag ratio for several vertical center-of-gravity locations. Stability derivatives used in 
the calculations were based on the results of static- and dynamic -force tests of a para- 
wing model having an aspect ratio of 2.83 and a deployed leading-edge sweep angle of 50°. 
The effects of center-of-gravity location on these stability derivatives are considered 
important in themselves to the understanding of the dynamic stability and control of para- 
wing vehicles. Consequently, the effects of center-of-gravity location on the stability 
derivatives are analyzed as a distinct part of the investigation. 

SYMBOLS AND NOMENCLATURE 

The calculated stability and control results are presented with respect to the sta- 
bility axis system shown in figure 1 . All basic force test data of the parawing model are 
presented for the bocty axis system shown in figure 2. The aerodynamic coefficients are 
based on the flat-pattern characteristics of the parawing (45° leading-edge sweep condi- 
tion) . Indicated dimensional units are given first in the International System (SI) and 
parenthetically in U.S. Customary System of Units. 

A,B,C,D,E coefficients of lateral stability quartic 

b wing span, meters (feet) 

c mean aerodynamic chord, meters (feet) 


2 


Cd 

Cl 

Cz 

AC^ 

Cm 

Cn 

AC„ 

Cy 

ACy 

Cl/2 

I>b 

f 

Fd 

Fl 

Fy 

g 

k 

^Xo 

^ZO 


drag coefficient, 
lift coefficient, FL/q„<,S 
rolling-moment coefficient, Mx/q«jSb 
incremental rolling- moment coefficient 
pitching-moment coefficient, 
yawing-moment coefficient, M^/^ooSb 
incremental yawing- moment coefficient 
side-force coefficient, Fy/qoo^ 
incremental side-force coefficient 

number of cycles required for amplitude of a lateral oscillation to 
decrease by a factor of 2 

differential operator, d/ds^, 

frequency of oscillation, cycles per second 

drag, newtons (pounds) 

lift, newtons (pounds) 

side force, newtons (pounds) 

acceleration due to gravity, meter s/second^ (feet/second^) 
reduced frequency parameter, ^ 

radius of gyration in roll about principal longitudinal axis, meters (feet) 
radius of gyration in yaw about principal vertical axis, meters (feet) 
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^Xo 




Xx 


K 


Z 


Xxz 

L/D 

m 

Mx 

My 

Mz 

P 

P 

^loo 

r 

S 


Sb 


nondimensional radius of gyration in roll about principal longitudinal axis, 



nondimensional radius of gyration in yaw about principal vertical axis, 

“Zo/” 

nondimensional radius of gyration in roll about longitudinal stability axis, 

IZ 2 2 2 . 2 

yXxo cos r] + Kzq sin 77 

nondimensional radius of gyration in yaw about vertical stability axis, 

^XxQ^sin2T7 + Kzq^cos^tj 

nondimensional product-of-inertia parameter, ” XxQ^jsin rj cos 77 

lift- drag ratio 
maximum lift- drag ratio 
mass, kg (slugs) 

rolling moment, meter-newton (foot-pound) 
pitching moment, meter- newton (foot-pound) 
yawing moment, meter-newton (foot-pound) 
period of oscillation, second 
rolling velocity, radians/second 

free -stream dynamic pressure, newton/square meter (pound/square foot) 

yawing velocity, radians/second 

wing area, square meter (square foot) 

nondimensional time parameter based on span, Vt/b 
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time, second 


time reqvdred for amplitude of lateral oscillation to decrease by a factor 
of 2, second 

free-stream velocity, meters/second (feet/second) 

^ « V 'fcp 

body reference axes 

nondimensional distances along the body reference axes measured from 
basic moment reference center (see fig. 3) 

nondimensional distances measured from 0.5 keel length station to center- 
of-gravity location measured in axis system parallel to stability axes 
(see fig. 1) 

any of mass or aerodynamic terms in lateral equations of motion 

angle of attack of keel, degrees 

angle of sideslip, degrees or radians 

rate of change of sideslip, radians/second 

angle of sideslip of parawing, degrees 

angle between reference axis and principal axis, positive when reference 
axis is above principal axis at nose, degree (see fig. 1) 

angle of attack of principal longitudinal axis, positive when principal axis 
is above flight path at nose, degree (see fig. 1) 

flight-path angle, degree (see fig. 1) 

root of stability quartic AX^ + BX^+ CX^ + DX + E = 0 


lateral relative- density factor, m/pSb 


<p angle of bank, degrees or radians 

(pfff angle of bank of parawing, positive when right wing tip is down, degree 

a ratio of air density at altitude to that at sea level 
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angle of yaw, degrees or radians 

ratio of amplitudes of roll and yaw present in a mode of motion 
ratio of amplitudes of roll and sideslip present in a mode of motion 


ratio of bank-angle amplitude to equivalent-side-velocity amplitude for 


oscillatory mode. 



p mass density of air, kilograms/cubic meter (slugs/cubic foot) 

w = 27Tf, radians /second 


aC; 


per degree or per radian 


aCr 


Cno = — per degree or per radian 

P 0|3 


9Cy 

Cyn = — degree or per radian 
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aC; 

^2V 


Cnv^ - 
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Cnf, - 


9 

4v2 


% 

% 


9Cy 

pb 

^Iv 


9Cy 

9^ 

4V" 


6 





rb 

^2V 


CYr 


dCy 

rb 

^2V 


C,. = 


9Ci 
_ rb" 


Cn;,= 


dC 


n 


rb" 


^Yi- = 


9Cy 
. rb^ 


4V" 


4V" 


4V" 



dCi 



c„ ■ = 

^13 


8C 


n 


2V 



9Cy 



Dots over symbols denote derivatives with respect to time. 

In the present investigation the term "in-phase derivative" refers to any one of the 
oscillatory derivatives that is based on the components of forces and moments in phase 
with the angular displacement during the oscillatory tests. The term "out-of-phase 
derivative" refers to any one of the stability derivatives that is based on the components 
of forces and moments 90° out of phase with angular displacement. The oscillatory 
derivatives of the present investigation were measured in the following combinations: 


Rolling oscillation tests: 


C7 sin a - k^Cy . 

^/3 ^p 

Cn sin a - k^C„ . > 

^/3 ”p [ 

Cv sin a - k^Cv 

^|3 ^p J 


In phase 


CZp + sin a 

Cnp + Cn^ sin a > Out of phase 

CYp + Cy 0 sin a 
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Yawing oscillation tests: 


C? cos a + k'^C/ . 

Cn^ cos a + k^Cni- ') phase 


Cy cos a + k'^Cy. 


C? -Cl • cos a 

Ij. 


j 


Cnr - <^n^ cos « 


/ Out of phase 


^Yr " ^Y^ cos a 


METHOD OF ANALYSIS 

The analysis was made for a hypothetical paraglider system consisting of a 
parawing- payload combination having the mass characteristics of a recently proposed 
recovery system. The parawing and payload were assumed to be rigidly joined so that 
there was no relative motion between them except for intentional control movements. 
Longitudinal trim was assumed to be obtained solely by varying the fore-and-aft position 
of the center of gravity of the system. The vertical center -of -gravity location was varied 
by locating the payload relative to the parawing in such a manner as to yield several 
specified center-of-gravity locations perpendicular to the parawing keel member (in 
terms of z/b) while maintaining the proper position parallel to the keel (in terms of x/b) 
required for longitudinal trim. The stability derivatives for each center-of-gravity loca- 
tion were then estimated from measured derivatives for the wing alone, and calculations 
were made to determine the dynamic lateral stability and control of the vehicle for each 
center-of-gravity location. 


Description of Vehicle 

The recovery system was assumed to employ an aspect-ratio- 2.83 parawing having 
a conical canopy and flat- pattern sweep angle of 45°. In the deployed or flight condition, 
the parawing fabric was supported in a 50° sweep condition by three rigid tubular mem- 
bers which formed the keel and leading edges of the parawing. Parawing and payload 
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weights were 1779 newtons (400 pounds) and 15 124 newtons (3400 pounds), respectively. 
Wing loading was 337 newtons per square meter (7.04 pounds per square foot). The pay- 
load aerodynamics were assumed to consist only of drag and were considered only in the 
determination of longitudinal trim. The system center- of- gravity location was varied in 
the previously discussed manner such that three different configurations of the original 
concept were studied. These configurations are herein referred to as A, B, and C 

= 0.25, 0.50, and 0,75, respectively^. 

The configurations were assumed to use the wing-bank system for lateral control; 
that is, the wing was banked about an axis parallel to the keel member to produce the 
forces and moments required for lateral control. This control system is, in effect, 
similar to the center -of -gravity shift type of control actually used on parawings with sus- 
pended payloads. As stated in reference 3, when such a system is used for control the 
incremental lateral force and moment coefficients produced by wing bank may approxi- 
mately be expressed as 


ACy = *^w 

ACn = Cn^^w + I Cl sin S (1) 

^C^ = sin 


where sin /3^ = sin a sin (p^. 

The first terms of the right-hand side of the equations arise from the fact that a 
value of sideslip is obtained at the parawing when the wing is banked for keel angles of 
attack other than zero. The second terms are the contributions of the lateral component 
of the lift vector which has been tilted by banking the wing. 

Determination of Stability Derivatives 

The lateral stability derivatives of the configurations were assumed to be due to the 
parawing alone. There was however virtually no information available on the dynamic 
stability derivatives; therefore, in order to provide reasonable inputs for the dynamic 
stability analysis, an experimental investigation was conducted to determine both the 
static and dynamic stability derivatives of a 0.15-scale model of the parawing assumed 
for the theoretical analysis. 
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A sketch of the model parawing is shown in figure 3 and the geometric character- 
istics are listed in table L The model was constructed of nonporous plastic membrane 
attached to three equal length rigid members such that in the deployed or flight condi- 
tion, the rigid members supported the fabric in a 50° sweep condition. Additional infor- 
mation pertaining to the model can be found in reference 5. 

The static- and dynamic-force tests were conducted in a low-speed tunnel with a 
12-foot octagonal test section at the Langley Research Center. The tests were made at a 
dynamic pressure of 78.1 newtons per square meter (1.63 pounds per square foot) which 
corresponds to a Reynolds number of 0.99 x 10® based on keel length. The static lateral 
tests were made for an angle of sideslip range of ±5°. Forced oscillation tests were 
made in yaw and roll with angular amplitudes of ±5° and for a value of the reduced fre- 
quency parameter k of 0.25. A detailed description of the dynamic-force test equip- 
ment and the method of obtaining the dynamic parameters is presented in reference 6. 

The moments were all measured about the basic moment center shown in figure 3; and 
the static and dynamic lateral stability derivatives of the complete configurations were 
obtained by transferring these measured wing-alone model data to the required center- 
of-gravity locations by means of the equations of reference 7 which are included in 
appendix A of this paper for convenience. 

Method and Scope of Calculations 

Calculations were made to determine the dynamic lateral stability and control of 
the three parawing configurations for the condition of (L/D)max gliding flight at an 
altitude of 3048 meters (10 000 feet). 

The lateral stability calculations consisted of the determination of the damping and 
period of the lateral modes of motion. The damping of both the oscillatory and aperiodic 

modes is expressed in terms of the damping factor — — , the reciprocal of the time to 

h/2 

damp to one-half amplitude. Positive values of this parameter indicate stability (posi- 
tive damping) and negative values indicate instability (negative damping). The calcula- 
tions were made with the use of the equations of motion given in appendix B. 

Calculations were also made to determine the aerodynamic and mass parameters 
that affected the lateral stability of the configurations to the greatest extent as the center- 
of-gravity location was varied. These calculations were made by the method of refer- 
ence 8 which is based on a Taylor' s series expansion of the roots of the stability quartic 
near the original solution. Although strictly applicable only for small incremental 
changes in derivatives and mass parameters near the original solution, the method does 
afford an insight into the important terms in the equations of motion as the center of 
gravity, and hence the stability derivatives and mass distribution are changed. The 
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results are presented as slopes indicating the rate of change of the damping and frequency 
of the various modes of the lateral motion with changes in any of the parameters in the 
equations of motion. In an attempt to provide a better understanding of the magnitude of 
the factors causing the effects of vertical center- of -gravity position shown by comparison 
of the results for configurations A, B, and C, additional calculations were made for con- 
figuration A in which the values of mass and aerodynamic parameters were changed one 
at a time to the values for configuration B to determine the effects of these individual 
parameter changes on damping and period. 

The lateral control calculations consisted of calculation of the lateral motions 
following a step input of 5° wing bank angle. These calculations were made with a digital 
computer and used the basic equations of motion of appendix B together with control 
inputs calculated from equations (1). 

RESULTS AND DISCUSSION 

The results of the parawing model force tests are presented in figures 4 to 7. 
Calculated lateral stability derivatives and mass characteristics of the three full-scale 
configurations are presented in figures 8 and 9 and are given in table n. The results of 
the dynamic lateral stability and control calculations are presented in figures 10 to 14, 
the sensitivity of the dynamic lateral stability of the configurations to changes in the sta- 
bility derivatives and mass parameters being given in tables in and IV. The two primary 
results of the investigation - the changes in the stability derivatives and mass param- 
eters with center-of-gravity location and the effects on dynamic lateral stability and con- 
trol of these parameter changes - are discussed separately. 

It should be noted that any theoretical analysis of the dynamics of a parawing con- 
figuration can be limited by the concept itself. The flexibility of the parawing, the means 
of connection between parawing and payload, and the number of degrees of freedom of the 
system all complicate and in some cases may render useless any analytic attempt to pre- 
dict the dynamics of the system. Consequently, the calculated results contained herein 
do not apply directly to parawing configurations or flight conditions other than those 
investigated. It is believed, however, that the major results and trends shown by these 
calculations are indicative of the dynamic characteristics to be expected of vehicles which 
have the following characteristics: very low center-of-gravity location, mass distributed 
mainly along the vertical axis, and a wing which has a large value of the lateral-force 
derivative Cy^. These characteristics are typical of a wide variety of paragliders or 

gliding parachutes. 
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stability Derivatives and Mass Parameters 

Mass distribution.- The calculated mass and geometric characteristics of the three 
configurations are given in table n. It should be noted that for all three configurations, 
the nondimensional radius of gyration in roll was larger than the nondimensional 

radius of gyration in yaw This situation is very different from that of a conven- 

tional airplane where the radius of gyration in roll is less than that in yaw < K 2 and 
leads to near-vertical inclinations of the principal axis of least inertia. This difference 
would be expected to lead to substantial differences in the character of the lateral motions 
between parawing vehicles and conventional aircraft. For example, the large values of 
inclination of the axis of least inertia will lead to increased yawing components in the 
lateral modes of motion since a vehicle will inherently tend to move about an axis of 
least inertia. 

It should be noted that in the present investigation the principal longitudinal axis 
(the axis to which 77 is referenced) was considered, in the conventional manner, to be 
the axis located in a generally forward direction. The only unusual situation was that for 
the parawing vehicles this axis was the one of greatest inertia. 

Longitudinal characteristics.- The para wing system was assumed (on the basis of 
the model force-test data of fig. 4 with an added drag increment to account for the pay- 
load) to have a maximum L/D value of 4.27 at a keel angle of attack of 25° for all 
center -of- gravity locations. The longitudinal stability of the configurations was not 
investigated in this study, but it might be noted that a paraglider configuration becomes 
more stable longitudinally for increasing values of x/b and z/b. 

Static lateral stability derivatives,- The static lateral stability characteristics of 
the parawing model are presented in figure 5 in the form of the static stability deriva- 
tives Cy^, plotted against angle of attack. The values of the derivatives 

were obtained from the differences between the values of the coefficients at sideslip 
angles of 5° and -5°. The variation of (referenced to body axes) shows a negative, 

or destabilizing trend at an angle of attack of about 33°. It should not be inferred from 
this result that a vehicle with this parawing would experience a directional divergence at 
these higher angles of attack. It is a well-established fact that when an aircraft has large 
negative values of in conjunction with the negative values of a directional 

divergence does not necessarily occur. This point is discussed in detail in reference 9. 

The values of the static lateral stability derivatives (referenced to stability axes) 
for the three configurations based on the model data and the transfer equations of the 
appendix are given in figure 8 (a). Lowering the center -of -gravity location from configu- 
ration A to that of configuration C increased the values of both and The 
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increase in the directional stability derivative Cn^ is caused by the fact that as vertical 

center- of- gravity location zA> is increased, the horizontal location x/b must also be 
changed to maintain longitudinal trim. 

Dynamic lateral stabili ty derivatives .- The variation of the lateral in-phase oscil- 
latory derivatives with angle of attack for the parawing model is presented in figure 6. 
Static data referenced to body axes are also presented and are shown to be generally in 
good agreement with the oscillatory data. 


The variation of the out-of-phase derivatives with angle of attack is presented in 

figure 7. The damping-in-roll parameter +C^*sino! and the damping- in-yaw 

P p 

parameter cos a were essentially constant over the angle- of -attack range up 

to the stall, which probably indicates that the ^ portion of the derivatives were negli- 
gible in the angle-of-attack range below about 35°. This apparent lack of ^ derivatives 
is in contrast with results of forced oscillation tests of swept rigid wings (see ref. 10, 
for example) in which separated flow effects have led to large values of the /3 deriva- 
tives at high angles of attack. The contrasting characteristics of the parawing may have 
been caused by the flexibility of the parawing itself or by leading-edge characteristics. 

In any event, it was assumed that the rate of change of sideslip (|3) derivatives were neg- 
ligible and the measured combinations were treated as pure rolling and yawing deriva- 
tives in the dynamic stability calculations. For example, the combination C^p + sin a 

was assumed to be for purposes of determining the lateral stability derivatives of 

the full-scale configuration. 


The calculated effect of center-of-gravity location on the dynamic lateral stability 
derivatives of the three full-scale parawing configurations is shown in figures 8(b) 
and 8(c). Lowering the center-of-gravity location from that of configuration A to con- 
figuration C led to a positive increase in the yawing moment due to rolling velocity Cn„ 

and negative increases in the damping-in-roU derivative C^p and the damping-in-yaw 
derivative The changes in the dynamic stability derivatives with vertical center- 


of-gravity location are very large because of the combination of large distance between 
the center of gravity and the wing and large lateral-force derivative Cy^ which enter 

into the transfer of these derivatives as products. In effect, increasing the vertical dis- 
tance between the center of gravity and the parawing keel essentially makes the parawing 
act as a high vertical tail (with short tail length) would on a conventional airplane — that 
is, the side-force derivative Cy^ of the parawing tends to be the predominant contribu- 
tor to the derivative transfer equations. 


13 


Dynamic Lateral Stability and Control 


The results of the calculations to determine dynamic lateral stability are discussed 
in terms of the stability of the different modes of lateral motion. For conventional air- 
craft, these modes are usually a highly damped aperiodic rolling motion known as the 
roll-subsidence mode, a lightly damped aperiodic motion involving yawing and rolling 
known as the spiral mode, and a lateral oscillation involving yawing, rolling, and side- 
slipping known as the Dutch roll. Although the lateral modes of a parawing configuration 
may be slightly different in nature in comparison with conventional aircraft modes, the 
conventional nomenclature will be used to identify the various modes. 


Roll-subsidence mode.- The calculated damping factor of the roll- subsidence mode 
of the configurations is presented in figure 10. The data indicate that the system became 
less damped as the center -of -gravity location of the system was moved downward from 
that of configuration A to that of configuration C. The dashed line represents the classi- 
cal single-degree-of -freedom approximation to the roll-subsidence mode damping 

/ 1 




(tl/2 1.57MbKxV 


which generally gives a reasonably accurate approximation for 


conventional aircraft. Comparison of the two results indicates large differences between 
the complete three-degree-of-freedom results and the simple single-degree-of-freedom 
results. This result indicates that the damping of this mode of motion is dependent on 
factors other than those included in the simple approximation. 


The calculated variation in the damping factor of the roll- subsidence mode with 
changes in the various parameters in the equations of motion are presented in table HI. 

o 1 


t 


Positive values of 
of 


1/2 


dXi 


together with positive increases in or negative values 


1/2 


axi 


together with negative increases in xW lead to positive (stabilizing) increases 


in the damping of this particular mode of motion. The results show that the yawing- 
moment derivatives and inertia pal-ameters can appreciably affect the damping of this 
mode. The results of the additional calculations in which the parameters in the equa- 
tions of motion of configuration A were changed one at a time to those of configuration B 
in an effort to determine the relative magnitude of the changes in the roll-mode damping 
due to changes in mass and aerodynamic parameters are presented in table IV together 
with results predicted by the slope method. The results show that the decrease in 
damping of the roll mode as the center-of-gravity location changed from that of configu- 
ration A to configuration B was caused primarily by the destabilizing contributions of 
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the derivative and the mass parameter . The increased importance of the 

yawing derivatives in the determination of the damping of the roll mode can be explained 
by the data of figure 11 which presents the ratio of roll to yaw amplitude in the roll mode. 
Values of this parameter usually are in the range of 30 to 100 for conventional aircraft 
but for these parawing configurations with their associated stability parameters, all 
values were below 5 and indicate that the mode was not a pure rolling motion but involved 
relatively large yawing motions. In fact, for the lower center- of-gravity location (con- 
figuration C), the motion consisted of almost an equal rolling and yawing - which is very 
different from the case for conventional aircraft where the motion is almost pure rolling. 

Spiral mod e.- The calculated results presented in figure 10 indicate that the spiral 

mode became slightly more damped as the center of gravity was lowered with respect to 

the wing. An inspection of table IV indicates that the stabilizing effect was primarily due 

to the stabilizing influence of changes in those derivatives and that usually 

affect the spiral stability of an aircraft. The spiral root can usually be approximated by 
E 

X = - — and those terms which most affect this ratio are seen to be of great importance 
in the damping of the spiral mode for the parawing configuration — as they have tradi- 
tionally been for conventional airplanes. 


Oscillatory mode .- Presented in figure 12 is a summary of the Dutch roll oscilla- 


tion characteristics of the three configurations. The inverse cyclic damping 


Cl/2 


and 


roll-to-sideslip parameter |— 1 for each configuration are presented together with the 

|Ve I 

recommended boundary of reference 11 in figure 12(a). The presentation of the Dutch 
roll characteristics in this manner is not necessarily meant to imply that parawing con- 
figurations should conform to the handling- quality parameters of conventional aircraft, 
but it is one of the present standards available for evaluating the effect of the Dutch roll 
characteristics on flying qualities and is consequently used for purposes of comparison. 
The results are also compared with the older handling- quality requirements (expressed 
in terms of damping and period of the oscillation) of reference 12 inasmuch as the 
present-day aircraft requirements tend to present pessimistic views of the oscillatory 
characteristics of vehicles having extremely light wing loadings. (The equivalent side 
velocity Ve of such vehicles is much smaller than that of conventional aircraft.) The 
results indicate that increasing the vertical distance between parawing and payload 
(varying center- of- gravity location from configurations A to C) led to a longer period, 
better damped lateral oscillation, and smaller values of roll-to-sideslip ratio. Only con- 
figuration C satisfied the present-day handling- quality requirements, both configura- 
tions A and B lacking satisfactory Dutch roll damping. The results presented in fig- 
ure 12(b) show that only configuration A did not satisfy the older handling- quality 
requirements. 
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Tables in and IV indicate that the increase in damping as the center of gravity is 
lowered is primarily caused by the stabilizing contribution resulting from the positive 
increase in Cnp being larger than the destabilizing contribution of The effects 

of Cnp on the damping of the lateral oscillation and aperiodic modes have been inves- 
tigated in the past (refs. 13 and 14) where it was shown that positive increases in the 
derivative usually increases the damping of the Dutch roll mode while decreasing the 
damping of the roll- subsidence mode. The derivative Cnp will almost always be 
increased positively by increasing vertical distance between the parawing and center of 
gravity. Increases in the same geometric variables will usually lead to negative values 
of Kxz explained previously. 

It should be noted that in table III the contributions of Cn^j C^^, and C^^ to the 
Dutch roll damping reverse signs between configurations A and B; that is, increases in 
the derivatives in the usual stabilizing sense ^negative for C^p, positive for and 

Cl ^ lead to less stable Dutch roll oscillations for configurations B and C. This result 
is due to the effect of C^p on the distribution of damping of the system. Similar 
results are reported in reference 13 for large positive values of Cnp. As discussed in 
reference 14, the total damping of the system can be expressed as B/A. The partial 
derivative of the ratio with respect to any parameter xj must therefore follow the 
relation 


0.693b 9xj y 9xi /Spiral 

mode 


9xi /Roll 

mode 


1 /Dutch roll 
mode 


Since for Cn« and C; ^ must be zero (changes in these derivatives do not 

9xj P 

change the total damping but merely redistribute the damping among the modes), the 
stabilizing contribution of the derivatives to the roll- subsidence damping is accompanied 
by a decrease in the Dutch roll damping. 

Theoretical investigations in the past (see ref. 1 5, for example) have shown that the 

L 1 \ 


reversal of the sign of 


\^^/iyntch roll 
mode 

dieted by the algebraic sign of the expression 


between configurations A and B could be pre- 


2ClKz^ 
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Positive values of this quantity indicated that increases in the effective dihedral param- 
eter -Ci^ would lead to larger values of the Dutch roll damping factor. The reversal 
of the sense of the contribution of ^ to the damping of the lateral oscillation as the 
center of gravity was lowered with respect to the parawing was therefore due to the posi- 
tive increase in Cnp brought about by center-of-gravity position. 

The data presented in tables HI and IV indicate that the decrease in frequency of 

the Dutch roll oscillation with increasing values of x/b was due primarily to the mass 

2 

parameters Kxz and Kx and the stability derivatives Cn^, Cnp, and 

Although changes in lateral stability with center-of-gravity location have been dis- 
cussed only for changes from that of configuration A to that of configuration B, similar 
results occur in going from configuration B to configuration C. 


Lateral response to wing bank control .- The incremental lateral-force and moment 
coefficients presented in figure 9 indicate that as the center-of-gravity location was 
changed from that of configuration A to that of configuration C, the incremental rolling- 
moment coefficient AC^ increased as a result of the increased value of z/b. The 
incremental yawing- moment coefficient ACn is negative, or adverse, for all configura- 
tions. The responses of the configuration to 5° wing bank control are presented in fig- 
ure 13. The results of the three-degree-of-freedom calculations show that the initial 
roll response was about equal for all configurations, but as the motion progressed, the 
configurations with the lower center-of-gravity locations (configurations B and C) showed 
greater response to wing bank control. The enlarged plot of the time history of the ini- 
tial yawing motions shown in figure 14 indicates that increasing the vertical center-of- 
gravity location decreased the initial adverse yawing motion, which is opposite to what 
would be expected from a static aerodynamic standpoint. This result is caused by the 
favorable contribution of the product-of-inertia factor to the initial yawing acceleration 
as expressed by the relation: 



ACn 


^XZ 

Kx2 


AC^ 







( 2 ) 


As pointed out in reference 16 and as can be seen from equation (2), negative values of 
*^xz give a favorable yawing-moment contribution and can thereby significantly affect 
the yawing characteristics in an aileron roll as the center-of-gravity location is lowered 
from that of configuration A to that of configuration C. The product-of-inertia parameter, 
as discussed earlier, will almost certainly be negative for large values of z/b. 
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SUMMARY OF RESULTS 


The results of a theoretical investigation of the dynamic lateral stability and con- 
trol of a parawing configuration having rigid leading e<%e and keel members may be smn- 
marized as follows: 

1. Increasing the vertical distance between the center of gravity and the parawing 
keel while maintaining the horizontal location required for longitudinal trim led to a 
decrease in the damping of the roll-subsidence mode and an increase in the damping of 
the spiral mode. Lowering the center of gravity also caused an increase in the damping 
and period of the lateral oscillation and a reduction in the ratio of roll to sideslip. 

2. The two parameters that most affected the damping of the lateral oscillation 
with vertical changes in center-of-gravity location were the yawing-moment-due-to 
rolling derivative Cnp and the nondimensional product-of-inertia factor 

3. There was little effect of vertical center-of-gravity location on the initial roll 
response to wing bank control but as the motion progressed, the configurations having 
the lower centers of gravity displayed the greater roll response. The effects of vertical 
center-of-gravity location on the mass distribution of the system produced significant 
changes in the initial adverse yawing motion following the application of wing bank control. 
Negative values of Kxz caused by lowering the center of gravity relative to the para- 
wing tended to reduce the initial adverse yaw. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., February 24, 1966. 
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APPENDIX A 


STABILITY DERIVATIVE TRANSFER EQUATIONS 


The stability derivative transfer equations given herein were taken from refer- 
ence 7 and are expressed in the notation of this paper. The distances x and z cor- 
respond to those of figure 3, Zero subscripts indicate data measured about the basic 
reference center. The equations convert body-axis derivatives to stability-axis 
derivatives. 

Cyo = Cv 




Cy„ = cos O' + C„ sin O! + 20^. 

^pq yro yj3o\l> 


cos 0 ! - ^ sin a 
b 


Cv = Cv cos a - Cv., sin o- - 2Cvo cos o; + § sin O') 
^r ^rQ ^PO ^Po\h b } 


Cj Cj ^ cos O' + Cn ^ sin a + Cv„ cos a - ^ sin 

^Po Po t) 


a 


C? = C; cos^Q!+ fc? + Cn„ \sin O' COS Q!+ Cji sin a 
^P ^Pq ( ^ro ^Pq) "^^0 


Cv + 2C? \cos a + (Cv + 2Cn„ \sin a\ 

% ^PO I ^^0 %' * 


— cos O' - — sin O' 


+ 2Cv (— cos O' - ^ sin a 


Cj = C; COS^O! + (Cn^ - C; \sin a cos a - C« sin^o? 
^r Itq ^ % ipoj “Po 


2C^ cos^a + [Cy + 2C„ ]sin a cos a - Cy sin^o; 


+ 2Cy 


Po “^o; "Po 

Cv COS^O! - fcy., + 2C; ^ \sin o; cos a - 2Cno sin^o? 

^PO ^0 

:j 


Po 


\b2 b2 


sin o! cos a - —icos^a - sin2o!j 


III ill III Hill III I III II III! Ill II 1 1 II II mil I III I III nil 


II mil Miiiiiiiiiiiiiiiii 
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Cno = Cno COS a - Cl ^ sin a - Cv^ Ir cos or + ~ sin o') 
p Pq ‘•/So b / 

Cn„ = Cnn_ cos2q! + J^sin a cos a - C7 sin^o- 


P l^-“ro -‘Po/ 

Cv cos^Q! + (Cv + 2 Cno Vin a cos a - 2C; sin^a 

^Po I ^0 ^oj ^^0 


2Cn^ cos^Q! - (Cv + 2C? \ sin O' cos O' - Cv sin^O! 

[ ^Po ^0 


+ 2Cy 


% 


K 

z2\ . 

:r Sin a cos a - 

^(cos2q! - sin2o!) 

Vb^ 

b2/ 



Cn^ = Cji^ cos^Q! - + Cn^ ^sin a cos q!+ C^.. sin^a 


%)=*■ 


PO 


•Yro " “ ■ (%o ^ “ 


Cv + 2 C„ \ 


+ ^ cos a + 
b 


/X Z 

+ 2Cv I — cos a — sin o') 
^ YaA u b ) 


^ 0 \ b 
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EQUATIONS OF MOTION 

The nondimensional lateral equations of motion (ref. 13), referred to a stability 
axis system (fig. 1), are: 

Roll: 

2>‘b(%V‘f * %zDbV) = * i . i 

Yaw: 

+ | CnpD^^0 + | Cnj.D^j<// 

Sideslip: 

2Mb {Ph^ ^b*^) ^ + I + Cl(^) + | CYj.D^jI// + (Cl tan y)i// 

ASv, 

When (PqB ° is substituted for 4), \}/Qe for i//, and iSqO for ^ in the 
equations written in determinant form, X must be a root of the stability equation 
(neglecting the zero or heading root) 

AX^ + BX^ + CX^ + DX + E = 0 (Bl) 

where 

A . 8Mb^(Kx^K22 . Kxz^) 

B = -2bi,2(2Kx2Kz2cY^ . Kx^C„^ . - K^z<^‘r ' Kxz'="p) 

+ Ky^Cj^Cy^ + 1 C„^C;p - K^z<^i^Cy^ - 

- CnpKxzCY^ - i C„pC;^ + KxzCn^CYp - K^^CYpC;^ - Kx^CYpCp^ -b %zCYpC, J 
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APPENDIX B 


D = - i Cn^C^pCy^ - + i CnpC; + MpCnpC^jg + 

- ^ - Zii^Ky^Cfi^Cj^ tan y + tan y + ^ ^tp^n^*2y^ 

- J CnpCtgCYr - i Cj^CngCyp + i C„^CjgCyp 

E ' I CL(Cn,C,^ - C;^C„^) y 1 Cl tan y(c,pC„^ - CnpC;^) 

The damping and period of a mode of motion in seconds are given, respectively, by 

the equations t, /« = - — and P = — — where c and d are the real and imagi- 

1/ z c V d V 

nary parts of the root of the stability equation. 
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TABLE I.- GEOMETRIC CHARACTERISTICS OF THE MODEL 


Aspect ratio: 

Flat 

Deployed 

Area: 

Flat 

Deployed 

Mean aerodynamic chord: 
Flat and deployed . . . 

Span: 

Flat 

Deployed 


2.83 

2.57 


1.14 sq m (12.27 sq ft) 
1.04 sq m (11.16 sq ft) 

. . 84.6 cm (33.33 in.) 

. . 179 cm (70.71 in.) 

. . 163 cm (64.28 in.) 


Root chord 

Sweep angle: 
Flat . . . 
Deployed 


127 cm (50.00 in.) 

.... 45.00 deg 
.... 50.00 deg 
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TABLE II.- MASS AND AERODYNAMIC PARAMETERS OF THE CONFIGURATIONS 
(ai 1 values are presented with respect to stability axe^ 


Parameter 

Configuration ] 

A 

B 

c 1 

^ 



-0.2556 

-0.2556 

-0.2556 

0.072165 

0.082705 

0.09325 



% 

-0.14287 

-0.208458 

-0.27405 

0.0701 

-0.06107 

-0.1923 

CYr 

0.00314 

0.02422 

0.04529 


-0.08946 

-0.17846 

-0.33479 

C;" 

0.03054 

0.04852 

0.07725 

^np 

0.009 

0.0486 

0.09898 



-0.01757 

-0.02452 

-0.0332 



0.01815 

0.04118 

0.079694 

Kz^ 

0.0163 

0.01697 

0.018045 

%Z 

0.000681 

-0.003266 

-0.00972 



3.183 

3.183 

3.183 

Cl 

0.757 

0.757 

0.757 

tan y 

-0.2338 

-0.2338 

-0.2338 

V 

1 

00 

o 

h-^ 

CO 

7°36’ 

9°12’ 

x/b 

-0.0636 

-0.1348 

-0.2060 

z/b 

0.25 

0.50 

0.75 
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TABLE m.- VARIATION OF THE LATERAL STABILITY OF THE THREE CONFIGURATIONS 
WITH CHANGES IN THE STABILITY PARAMETERS 


c^p. 

^np - 

CYp. 
Olr . 
Cnr - 
CYr 


Cy, 




*'xz 

Ctp. 

C„p . 

CYp. 
Cjr • 

Cnj. . 

CYr • 


Cn, 


Cy. 

y 

K, 




^XZ- 


Cnp 

CYp 

CZr 

Cnr 

CYr 

Cl/3 

Cn|3 

cy^ 

y • 


'^xz* 




Configuration A 


Configuration B 


Configuration C 


Spiral 

Roil 

subsidence 

Dutch roll 

a 

a-^ 

8^ 

au) 

h/2 

*1/2 

*1/2 

axi 

dXi 

8X1 

axj 


m A 




-0.24 

-14.059 

-0.97 

-0.13 

-0.47 

-20.82 

11.012 

-4.88 

0.00021 

-0.18 

0.091 

0.11 

-3.77 

3.36 

0.55 

0.20 

-7.45 

4.97 

-7.81 

1.46 

0.0034 

0.043 

-0.023 

-0.079 

-0.67 

-1.23 

0.95 

-0.83 

-1.32 

-1.82 

1.56 

13.88 

0.00059 

-0.016 

-0.29 

0.0034 

-0.25 

-0.089 

0.17 

-0.017 

-0.062 

-85.36 

1.87 

-7.63 

-1.92 

30.18 

-24.81 

-60.69 

-1.091 

-106.045 

77.57 

89.41 

n B 




-0.029 

-9.22 

0.10 

-0.11 

-0.074 

-21.34 

10.011 

-2.56 

0.000023 

-0.099 

0.050 

0.089 

-4.26 

4.55 

-0.84 

0.013 

-10.74 

10.52 

-8.70 

1.32 

0.0034 

0.049 

-0.026 

-0.078 

-0.95 

2.052 

-0.55 

1.51 

-2.38 

4.65 

-1.13 

16.35 

0.00076 

0.022 

-0.31 

-0.0042 

-0.17 

-0.50 

0.33 

-0.044 

-0.010 

-33.83 

2.54 

6.65 

-3.69 

38.61 

-26.91 

-57.82 

-1.49 

-61.62 

44.78 

59.82 

)n C 




0.29 

-6.57 

1.12 

0.096 

0.86 

-18,88 

7.95 

-0.86 

-0.00017 

-0.046 

0.023 

0.066 

-5.36 

5.93 

-1.36 

-0.16 

-15.74 

17.064 

-9.39 

0.73 

0.0026 

0.042 

-0.022 

-0,079 

-1.77 

5.15 

-1.69 

2.32 

-5.13 

14.48 

-4.68 

17.50 

0.00092 

0.035 

-0.35 

-0.0091 

0.019 

-1.045 

0.51 

-0.027 

0.13 

-14.43 

0.52 

7.45 

-6.89 

37.49 

-23.94 

-58.095 

-1.97 

-28.45 

23.58 

39.97 
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TABLE IV.- VARIATION OF THE LATERAL STABILITY OF CONFIGURATION A 
WITH CHANGES IN THE STABILITY PARAMETERS* 



K 2 

%z 


Cnp 


% 

Cnp 

CYr 


Spiral 

Roll 

subsidence 

Dutch roll 

= 

*1/2 

*1/2 

Aw 

A 

*1/2 

).04118 

0.01630 

0.00068 

-0.14287 

0.07216 

0.07010 - 

■0.08946 1 

0.00900 0.00314 

0.03054 -0.01757 

-0.0014 

-0.93 

-0.081 

0.059 

■ 











(-.0014) 

(-1.97) 

(-.18) 

(.043) 

1 

.01815 

.01697 

.00068 

-.14287 

.07216 

.07010 

-.08946 

.00900 

.00314 

.03054 -.01757 

-.0013 

.020 

-.039 

-.016 

1 











(-.0013) 

(.020) 

(-0.41) 

(-.017) 

= 

.01815 

.01630 

-.00327 

-.14287 

.07216 

.07010 

-.08946 

.00900 

.00314 

.03054 -.01757 

.0045 

.51 

-.30 

-.31 

= 

■ 

s 











(.0043) 

(.42) 

(-.35) 

(-.31) 

s 

.01815 

.01630 

.00068 

-.20846 

.07216 

.07010 

-.08946 

.00900 

.00314 

.03054 -.01757 

.040 

.078 

.054 

-.059 

" 











(.044) 

(.081) 

(.055) 

(-.062) 


.01815 

.01630 

.00068 

-.14287 

.08271 

.07010 

-.08946 

.00900 

.00314 

.03054 -.01757 

-.013 

-.018 

.14 

.015 

= 











‘ (-.014) 

(-.019) 

(.15) 

(.016) 


.01815 

.01630 

.00068 

-.14287 

.07216 

-.06107 

-.08946 

.00900 

.00314 

.03054 -.01757 

-.000037 

.024 

-.015 

-.012 












(-.000028) 

(.024) 

(-.015) 

(-.012) 


.01815 

.01630 

.00068 

-.14287 

.07216 

.07010 

-.17846 

.00900 

.00314 

.03054 -.01757 

.012 

1.29 

-.0046 

.073 












(.021) 

(1.25) 

' (.011) 

(.086) 


.01815 

.01630 

.00068 

-.14287 

.07216 

.07010 

-.08946 

.04860 

.00314 

.03054 -.01757 

-.051 

-1.02 

-.15 

.55 

s 











(-.019) 

(-.83) 

(-.19) 

(.44) 


.01815 

.01630 

.00068 

-.14287 

.07216 

.07010 

-.08946 

.00900 

.02422 

.03054 -.01757 

.000062 

.00072 

-.0017 

-.00050 












(.000070) 

(.00092) 

(-.0017) (-.00049) 


.01815 

.01630 

.00068 

-.14287 

.07216 

.07010 

-.08946 

.00900 

.00314 

.04852 -.01757 

-.065 

.058 

.0036 

.0097 












(-.068) 

(.060) 

(.0036) 

(.0098) 


.01815 

.01630 

.00068 

-.14287 

.07216 

.07010 

-.08946 

.00900 

.00314 

.03054 -.02452 

.053 

-.037 

-.011 

.054 

I 











, (.051) 

' (-.035) 

(-.010) 

(.054) 

! 


^Derivatives changed one at a time to the values for configuration B. Derivatives chained are underlined. Variations predicted by 
the slope method of reference 8 are indicated in parentheses. 


! 



Reference axis (parallel to keel) 




Figure 1.-* The stability system of axes. Arrows indicate positive directions of moments, forces, and angles. This system of axes is defined as an 
orthogonal system having the origin at the center of gravity and in which the Z-axis is In the plane of symmetry and perpendicular to the relative 
wind, the X-axis is in the plane of symmetry and perpendicular to the Z-axis, and the Y-axis is perpendicular to the plane of symmetry. At a 
constant angle of attack, these axes are fixed in the configuration. 
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a, deg 


Figure 5.- Variation of the static lateral stability derivatives of the 0.15-scale model with angle of attack. Data referenced to body axes. 
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B 


C 


Configuration 


Static lateral derivatives. 
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